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ABSTRACT 


Turbulence data obtained over natural water waves were analyzed using 
joint probability distribution and conditionai means methods. These data 
represented conditions when the waves were decaying and when the waves 
were building. In both cases, significant wave-induced fluctuations were 
identified in the airflow. All features of the velocity fluctuations 
were examined for two levels above mean water level. In the case of a 
decaying wave field, decelerations in the airflow can be associated with 
an assumed propagating pressure maximum over the crest of the wave. Other 
than this deceleration, the airflow appears to reflect simple streamline 
bending over the mobile irregular wave surface. In the case of a building 
wave field, velocity fluctuations appear to agree with those predicted by 


linear wave generating theories. 
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I. INTRODUCTION 


This study examines turbulence data obtained over natural water waves 
in order to determine the influence of the waves on the overlying air. 

The data examined are simultaneous measurements of velocity and temperature 
fluctuations at two fixed levels above mean water level. 

The study extends a previous analysis on data for which significant 
evidence of wave-induced motion in the airflow was obtained. Joint proba- 
bility density and conditional mean analysis are the statistical procedures 
used to extend the previous interpretations. 

The lack of good descriptions of the interaction between the ocean 
surface and the adjacent air layer is a major obstacle preventing solution 
of a wide class of applied geophysical problems such as weather prediction 
and general circulation studies. A suitable description requires the 
application, to observational data, of specific statistical methods which 
could reveal properties of wave-induced organized motion. 

Current problems requiring additional knowledge on the dynamics of the 
surface layer are the determination of: 

1. kinetic energy dissipated by boundary layer turbulence; 

2. thermal energy gained by turbulent transport of sensible and 
latent heat; 

3. energy transformed by vertical velocities due to the horizontal 
distribution (convergence) of surface friction; and 

4. surface winds used in semi-empirical wave forecasting models. 

Although numerous observational investigations have been conducted to 
describe the structure of the near surface layer, the majority of these 


measurements have been over land. Therefore, parameterization of the 
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boundary layer over the sea has been based on overland results. It is 
clear, however, that a mobile irregular wave surface receiving kinetic 
energy from the airflow can influence the turbulent regime within the 
adjacent airflow. Therefore, empirical formulae, valid or at least 
tested over land, may not be valid over water. 

Examples of the conflicting results obtained from strictly empirical 
approaches appear in Figure 1 in which various proposed drag coefficient 
relations are presented for over the water. Being able to add som 
dynamical description to this region should help in narrowing ones' choices 
for this empirical coefficient. 

Data used in this study were previously analyzed by Davidson (1970). 
However, those analyses consisted of spectral analyses and the limitations 
imposed on the interpretation of spectral results hindered the description 
of several important properties of the observed crganized motion. Prece- 
dures employed in the present study are well suited for examining organized 
motion within a turbulent regime. The statistical procedures (Holland, 
1968) were designed to reveal systematic structure in various types of 
meteorological data but at the present time have only been applied to 
larger scale fluctuations. This study is the first attempt to apply these 


procedures to small scale fluctuations. 
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FIGURE l. 
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IT. THEORETICAL CONSIDERATIONS 


A. WIND WAVE-COUPLING 
In a survey of the understanding of the airflow over water, Stewart 

(1967) forewarned that the air-water boundary is not "just another 

boundary layer", but is unique in that the interface is a part of the 

solution for the boundary layer airflow. The interface is free to move 
and, as observed, it takes advantage of this freedom. Stewart, evaluating 
existing wave generation models (Miles and Phillips), indicated that in 

1967 there were still insufficient overwater measurements to either verify 

or invalidate existing analytical theories of the interaction between a 

turbulent shear flow and the underlying water surface. He reported, how- 

ever, that even though he examined a very large number of overwater 
velocity spectra, in no case were there sharply defined spectral peaks 
corresponding to the wave spectra peaks. Other studies, however, have 
revealed the existence of wave-induced motion in the air. 
Yefimov and Pososhkov (1969) reported: 
"Laboratory experiments indicate that surface waves induce in 
initially quiet air not only periodic velocity fluctuations put also 
stationary flow in the direction of wave propagation. The velocity 
of stationary flow is considerably higher than would follow if we 
assumed it to resemble Stokes wave flow in the water." 

Volkov (1969) reported, based on interpretations of spectral results: 
",.. the wave surface has a considerable influence on all of the 
statistical characteristics of turbulence in the atmospheric boundary 
layer - the fluctuation spectra and dispersion, the turbulent fluxes, 
etc. We therefore succeed in finding same typical features and 
regularities of this influence by the complex method of investigating 


the processes of interaction between the atmospheric boundary layer 
and the agitated sea surface." 


Tigh 





Davidson (1970), also interpreting spectral results, concluded: 


»-- analyses of the boundary layer over a majority of the earth's 
surface cannot be based on the assumption that the surface is a rigid 
boundary. The influence of the waves on the adjacent airflow has been 
found to be too significant to rely on over-land data to describe the 
boundary layer over waves. Reliable empirical relations giving the 
flux of momentum as a function of the wind profile and density gradi- 
ents will have to be formulated from data obtained over the ocean 
surface itself." 

An analytical model describing coupling between a wave surface and 
the overlying air flow was first proposed by Miles (1957). In this model, 
the wave-induced perturbation is considered in a coordinate system which 
is stationary with respect to the wave (i.e., surf-rider's coordinate 
system moving at the phase speed c). Streamlines in such a coordinate 
system would appear as in Figure 2. The Miles' critical height (dashed 
line) occurs at the level where the wind speed is equal to the wave speed. 
The wave, Figure 2, is stationary in this coordinate system; therefore, 
below the critical level the air appears to move to the ieft since it 
moves slower than the waves. Above the critical level, where the wind 
speed is greater than the phase speed, the air moves to the right. 

The Miles mechanism can be described as a positive feedback mechanism 
and, as such, it depends on the pre-existence of the wave. Forces within 
the Miles theory leading to wave growth are depicted in Figure 2. 

Region A is a region of negative vortex acceleration and Region B is one 

of positive vortex acceleration, but the absolute value of the vortex 

acceleration is greater at point A than at B. Therefore, there is a 

negative vortex acceleration at that level hence a deceleration of the 

mean wind when averaged over a wave length. This deceleration results in 


a loss of momentum from the airflow which contributes to the momentum 


increase of the wave. 
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FIGURE 2. Streamlines of sinusoidal perturbation in a 
shear flow near Z = Z_; (a) kinematic conse- 


quence, (b) consequente of dynamics of critical 
level. 
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It has been shown (Kinsman, 1965, p. 571) that the kinematics at the 
critical level, described above, would result in a phase shift, with 
height, of the wave-induced vertical velocity component. 

Therefore, it is observed, within a linear model, that the wave- 
induced fluctuations can be important in altering momentum transfer in 
the air layer adjacent to the waves. Also, specific predictions of this 
theory can be examined, such as a phase shift of the vertical velocity 
with respect to height. It should be noted, however, that within this 
investigation possible non-linear interactions between the wave-induced 


motion and the shear flow will also be examined. 


B. STATISTICAL CONSIDERATIONS AND PROCEDURES 

Statistical procedures used in this study were described by Holland 
(1968). Applied to three variables, these procedures yield the joint 
probability density function (JDF) for a pair of the variables with a 
conditional mean function (CMF) of a third variable. 

Holland argued convincingly on the suitability of this type of 
analysis for turbulent regimes. In reference to the application of these 
procedures to a turbulent regime, Holland cited Batchelor (1953): 

"It is a premise of probability theory that a random function f(a), 
say, defined for all values of a, is determined statistically by the 
complete system of joint-probability distributions of values of the 
function at any n values of a, where n may take any integral value. 
Likewise, the infinite field of turbulent motion is determined 
statistically by the complete system of jointzprobability distribu- 
tions of the values of the vector velocity u(x,t) at any n points of 
of space-time." 

A joint probability density function is a means to represent simul- 
taneously observed values of pairs of variables by contours on a two- 
dimensional array. A trivariate statistical relationship is obtained by 


computing a conditional mean of a third variable as a finction of the two 
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independent variables contributing to joint probability density function 
(JDF). This trivariate statistical relationship, hereafter denoted 
JDF/CMF, can therefore, be represented by two sets of contours on a 
two-dimensional array. 

Joint probability density function and conditional mean function 
analysis results are illustrated in Figure 3. The array is centered on 
the mean value, which is zero, of the two variables. Each cell in the 
row or columm consists of two numbers. The bottom number represents the 
JDF, the probability (x10?) of the joint occurrence of the representative 
deviation for the pair of variables, and the upper number represents the 
CMF, the mean of the third variable corresponding to the given joint 
occurrence of the pair. 

These results have been smoothed by a nine-point smoother to remove 
the perturbations viewed as being physically and statistically insignificant. 
This smoothing reduced an 18 by 18 array to the observed 16 by 16 array. 

The following list of steps within the JDF/CMF statistical procedures 
will be described in the ensuing paragraphs: 

1. normalization and computation of preliminary statistics of the 
variables; 

2. bandpass filtering; 

3. computation of JDF/CMF; 

4. checks for filter and random noise; and 

5. computation of phase-amplitude information. 

Computation and plotting of all statistics were done in the W. R. Church 


Computer Facility! of the Naval Postgraguate School. 


Ithe analysis utilized the IBM 360 Computer and a CALCOMP plotting 
device for the computations and contour plotting. 


15 





tO 
oOoo° 
ee @ 
ath 


uo a 
~ 


NO 
Ooo 

e@ee 
ich 


O 

Ls tt te tt 

2 Se Sx 4 

“WLU 

> a A 
AMM 


WO 


be 1 LL LU 
Ws Ld 
Mood 

qi 
<td pg 4 
m= Low 
<<aa¢ 
ao>>> 


Oo 


Ce 


oo oo 
e fe) 
Oo Lo) 

oo oo 
o oO 
© oO 

oO On 
e © 
oO © 
’ 

MAA AO 
@ e 
© oO 
i t 

MP NO 

e oN 
© © 


16 


Cutput Array 





R(XZ)=-0.08 RIYZ)}= 0.05 


R(XY)=-0.31 


= 63444 


COV(XZ)= -110.34% COV(YZ) 


COV(XY)= —-400.95 





1. Normalization and Computation of Statistics 

JDF/CMF statistics were computed from normalized values cf the 
Original variables. For example, if x (i) = values in one of the three 
sets, X = mean of all x(i), and o = standard deviation of all x(i), the 
JDF and CMF were computed for (x(i)-x)/o. Normalized values were obtained 
for each of the three variables included in a JDF/CMF analysis. 

In addition to the mean and standard deviation, other statistics 
were computed and tabulated for each variable and pairs of variables used 
in the analysis. These statistics were: 

l. variance; 

2. skewness of each variable; 

3. kurtosis of each variable; 

4. covariances; 

5. correlation coefficients: and 

6. number of data points used. 

2. Bandpass Filter 

A bandpass filter was applied in order to isolate the significant 
frequency band (0.1Hz to 0.8Hz) as indicated by previous spectral analysis 
on the data (Davidson, 1970). There were several numerical filters that 
could have been used. Investigation revealed that the most desirable 
filter, with respect to sharpness of cutoff and magnitude of oscillations 
beyond the terminal frequency, was an inverse transform filter. 

The response curve for the inverse transform bandpass filter with 
59 weights is shown in Figure 4. The filter is superimposed (Figure 14) 
on the wave spectra, also shown are the velocity variance and covariance 
spectra. The filter passes a range of frequencies from about 0.1 Hz to 


O23 Hz: 


17 





The application of this bandpass filter enabled analyses of wave- 
induced fluctuations without contamination by higher and lower frequency 


fluctuations. 


Bs ccisail 


| n (Hz) 1 10 


FIGURE 4. Response Curve for Inverse Transform Filter, 59 Pt. 


3. Computation of JDF/CMF 
A JDF represents the probability of a joint occurrence for a pair 

of variables. Therefore, for each pair of variables two indices were 
computed representing the position of simultaneously observed values of 
the variables in an 18 by 18 array of o/2 by o/2 joint class interval. 
For example, if the first values in the two sets were 1.25 and 0.9 
standard deviations from their respective means, their joint occurrence 
would be tabulated in the o/2 by o/2 cell with indices (3,2). The OF 
was computed for each joint class interval and was the average value of 
the third variable in the set for all observation times when the two 


independent variables had values within the joint interval. 


18 





An additional step in these procedures was to divide the number 
of joint occurrences by the total number of observations. This gave the 
probability per o“/4 within each cell of the 18 by 18 array which was 
then smoothed to a 16 by 16 array. Likewise, the sums of the third 
variable were divided by the number of joint occurrences in each cell. 
This yielded the average deviation from the mean of the third variable 
corresponding to a particular probability per o* for the other two 
variables. 

The output was printed out in a 16 by 16 array, Figure 3, and the 
array was analyzed by Contour, a computer library subroutine, and plotted 
on the CALCOMP plotting device. 

An analyzed JDF/CMF (for a random number field) is illustrated in 
Figure 5. The contours connect equal values of occurrences in the JDF. 
They are concentric circies because they were generated from a ranccm 
number field and circles indicate (as expected) no relationship between 
the two variables making up the JDF. A statistical relationship between 
two variables, other than between orthogonal relations such as a sine and 
cosine, would cause asymmetrical or elliptical contours. The heavier 
lines in this two-dimensional picture represent isolines of the CMF with 
the heaviest (dashed) being the zero line. An analysis with this degree 
of complexity, containing essential relations which exist among three 


variables, is not difficult to read with comprehension. 
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FIGURE 5. Analyzed JDF/CMF For Random Numbers 


Any pair of variables can be used to compute the JDF and any 
variable can be used to compute the CMF. In Figure 25, the JDF is u and 
w denoted by (u,w) and the CMF is n as a function of u andw. The JDF/ 
CMF is denoted by n (u,w). 

4. Checks for Filter and Random Noises 

Bandpass filtered data were analyzed and compared to the unfiltered 
data to insure that the data were not altered by the bandpass filter. The 
analysis of the same variables for the bandpassed (A) and the non-bandpassed 
(B) data are shown in Figure 6. It can be seen that the JDF (asymmetrical 
Circles) for the two sets of data are essentially the same. CMF features 
in the bandpassed analysis are smoothed versions of those in the non- 
bandpassed analysis. The resultant differences between A and B are differ- 
ences that could be caused by the addition of random numbers (which is 


essentially the contribution from background turbulence) onto the data. 
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FIGURE 6. JDF/CMF of Filtered (A) Data and Unfiltered (B) Data 


As a test on the results, random numbers were superimposed on the 
data used in each analysis. The values of these random numbers were of 
the same order of magnitude as the data field itself. It was expected 
that by superimposing random numbers, as in Figure 5, onto a data field, 
such as Figure 17A, where there is same statistical relationship, there 
would result a final JDF like Figure 17B where the asynmetry would be 
decreased. This was done for each variable analyzed and such a result 
was indeed realized. 

The analysis of each set of variables was done twice, once with 
the bandpassed data and once with a random number field superimposed on 
the bandpassed data. This allowed a comparison between the data and the 
data plus random numbers to examine if the significant modes within the 


results represented something other than chance occurrences. 
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5. Theory and Computation of Phase Amplitude 

In this study, information was also desired on the statistical 
dependence of the various turbulent variables on the phase of the waves. 
Within the context of a JDF/CMF analysis, this would be, for each variable 
plotted (Figure 18), its conditional mean function with respect to the 
phase angle of the wave. Similar to the JDF/CMF results a method of dis- 
play is required to emphasize particular variations and eliminate unwanted 
background noise. 

An approach suggested by Holland wes used in which the polar 
coordinates representing the amplitudes and phase angles of a variable 
can be determined by computing the JDF for a variable and its time deriva- 
tive. This method for computing phase relations should reveal (Holland, 
1968) if the phase relation between two variables shifts as a fimction of 
the amplitude of the reference variable. In particular, it will reveal if 
the amplitude of the dependent variable varies non-linearly with the ampli- 
tude of the independent variable. Such a result is highly significant in 
this study. 

Holland defined for any fluctuating meteorological variable y a 


"statistical eddy period", oe defined as 
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If the period of oscillation is t , the ratio of the standard deviations 


can be found. If the variable is written as 
y =A sin 2n= 
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The statistical eddy period for each variable, the standard devi- 
ation of each variable and the standard deviation of the time derivative 


of each variable are given in Table l. 










19 AUGUST : 26 SEPTEMBER | 
var. | 9, | | et | ea 


a503 | 8 [ase | an | an | 
uf 7687 [| 1002 | 336 | 39.92 | 140s | 3.4? 
W, 2.89 































Wu 
Tu | oos {oor [| 313 | oro | 0.04 | 3.36 _| 
= UW, 2.62 
UW, | 871.74 | 363.76 | 3.01 [1759.83 | 820.39 | 2.70 | 
wi 0.42 | 2.56 | 3.23 | 165 | 247 | 
wy | 0.87 [0.44 | 250 f 3.93 [83 | 2.70 | 


A 





TABLE 1. Statistical Eddy Period of Each Variable and Standard 
Deviations of Each Variable and its Time Derivative 
for 19 August and 26 September 1968 
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Phase results consisted of the statistical dependence of various 
turbulent variables on the phase of the waves (n). Therefore, a JDF 
coordinate system was used with n as the ordinate and its derivative (hn) 
as the abscissa. In this coordinate system, with the passage of a wave, 

n increases upward and successive values in time are higher on the right 
and lower on the left of the n-axis. 

If the surface wave were a pure sinusoid, successive values of 
JDF (n,n) would follow a counterclockwise circular path around the JDF 
coordinate system, passing through the extremes of n when crossing the 
N-axis and through extremes in fh when crossing the n-axis. However, water 
waves are not pure sinusoids so these paths would not necessarily be 
Circular. They would, however, be less complicated than passages of 
eddies in a random turbulent regine. 

The phase angie was defined py the angie measured (by octants) 
counterclockwise, Figure 7, from the positive n-axis in the (n,n) 
coordinate system. 

Phase relations were examined after separating the results into 
Wave amplitude classes. The classes were defined by a regions on the JDF 
array. In this case two amplitudes were considered, these being small 
amplitudes (less than 1.50) and large amplitudes (greater than 1.50). The 
numbered octants are shown in Figure 7 in which the box around the center 
represents 1.50. Essentially TS permitted an examination of phase 
relationship with respect to large and small waves. In contrast, spectral 
analysis yields phase information as a function of frequency only, the 


amplitude being the mean amplitude at that frequency. 
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SDF/CMF Smoothed Output Array Depicting the Octants fois 
Phase-amplitude Computations and 1.50 Boundaries (Box 
Around Center) for Separating Large and Small Amplitude 
Waves. 
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Applying concepts of the above discussion, computations of the 
phase amplitude relations were made using the dependent variables as the 
conditional mean and the waves and time derivative of the waves as the 
independent variables for the JDF. After obtaining the JDF/CMF smooth 
array output, for example Figure 3, the products of the dependent variable 
(the upper number in each cell) and the probability function (the lower 
number in each cell) in all joint intervals of each octant and amplitude 
class were computed to obtain the amplitude stratified phase conditional 
mean function. This was done, Figure 18 bottam to top, for n, u, w, T, 
-—uw, and wI, shown as horizontal bars for each octant. 

A smooth curve was drawn through each set of horizontal bars so 
as to make the two areas between each bar and the curve approximately 
equal. These smooth curves were then analyzed for their phase relation- 
ship with the wave curve (bottom curve). The left hand set of curves (A) 
represents the small amplitudes. The middle set (B) represents the large 


amplitude and the right hand set (C) is the total of A and B. 
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Iit. DESCRIPTION OF DATA 


The data examined in this study were from two days (19 August and 
26 September 1968) of measurements on Lake Michigan by Davidson (1970). 
These were simultaneous measurements from a fixed tower, Of wind (u and 
w) and temperature fluctuations and waves. The velocity and temperature 
measurements were obtained at two different levels above the mean water 


level. Waves were measured at a point directly below the air sensors. 


A. LOCATION AND DESCRIPTION OF THE TOWER? 

The measurements were made on Lake Michigan from a U.S. Lake Survey 
Research tower located in 15 meters of water 1.6 kilometers from shore 
near Muskegon, Michigan (Figure 8). Open water fetches at the tower 
location vary from &6 to 306 miles from the south clockwise through the 
northwest. Because the lake is not so large that separate weather 
systems and hence swell propagate from different directions, the local 
wave field was related to the wind field over the entire lake. 

The general construction of the research tower and its instrumentation 
is shown in Figure 9. The structure, an open triangular truss, was 
designed both to achieve stability and to minimize disturbance to the air- 
flow and water motions. These were important considerations in the present 
analysis and interpretation. Stability was achieved by supporting the 
tower on the lake bottom by a tripod which was secured by a concrete anchor 
at each leg. The tower was also guyed in order to withstand extreme wind 


and wave conditions. 


2an abridged version of the description by Davidson 1970). 
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B. INSTRUMENTATION AND SENSOR MOUNTING ARRANGEMENT 

Figure 10A shows the features of the sensor mounting arrangement. A 
plumb line arrangement of sensors enabled determination of phase relations 
between the waves and turbulent components in the air. 

The turbulent velocity measurements were made with a constant- 
temperature hot-wire anemometer using two wires in an X-configuration. 

The two wires arranged in X-configuration allowed the system to simul- 
taneously measure the two wind components of interest (u and w). The probes 
were calibrated both in situ, based on simultaneous cup and propeller 
anemameter measurements, and in a wind tunnel. 

A resistance thermometer, in a Wheatstone bridge circuit, was used to 
Measure the temperature fluctuations. The output signal of the system was 
125 millivoits for a one degree (Celcius) temperature change. Testing 
this system in a laboratory indicated it had the resolution to measure 
temperature to an accuracy of 0.05 C. 

The waves were measured by a capacitance gauge and bridge system, 
Figure 11. ‘The bridge was adjusted for full deflection, zero to one 
milliamp, between no immersion and full inmersion. 

The velocity and temperature sensors and wave gauge extended 
westward 2.5 times the width of the tower's side. Extension ams consis- 
ted of five centimeter rectangular tubing and were sturdy enough to avoid 
motion due to the airflow or waves. Also, a vertical pipe extended from 
below the water to two meters above the water so the sensors could be 
positioned at any height from the wave crest to two meters. Above the 
two meter level, arms were located at the four, eight, and fifteen meter 


levels. 
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Figure 10B is a photograph of the sensors mounted on the vertical 
pipe above the mean water level (MWL). A wind vane, such as that shown, 
was used at all levels to monitor the direction of the mean wind with 
respect to the sensors. The wave gauge can be seen directly below the 


hot wires in Figure 10B. 
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FIGURE ll. Picture of capacitance wave gauge and bridge- 
amplifier system. 
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IV. PRESENTATION AND INTERPRETATION OF RESULTS 


In this section, JDF/CMF and phase amplitude results are presented 
and interpreted for the most probable relations between fluctuations in 
the airflow and the waves. Temperature relationships were also computed 
and a cursory examination has revealed same interesting fluctuations in 
the temperature field, but a detailed analysis of the temperature relation- 
ships was not performed in this study. They are included for reference 
in Appendix B. 

Non-linear properties and phase relations of velocity fluctuations 
are emphasized in the following interpretation. Spectral results 
(Davidson, 1970) for these same data will be shown to complement the 
discussion. Spectral results are used to indicate the intensity of the 
observed wave-induced fluctuations with respect to the background, “normal” 
turbulent regime. 

Interpretations obtained from the analyses of the present study should 
be viewed as significant extensions to those available from the previous 
spectral results. Although the spectral results show, according to 
frequency, the existence of wave induced fluctuations, interpretations 
are limited with respect to non-linear properties of the fluctuations and 
phase relations between fluctuating variables. Non-linear properties are 
not available within spectral results because of the inherent linear 
expansion, i.e., Fourier Series. Also, phase information, obtained from 
spectral analysis, is highly uncertain because of the high level of 


background, random turbulence at all the frequencies. 
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For each period, two figures are used to show JDF/CMF results and 
two figures are used to show phase amplitude results. Each of the two 
figures representing JDF/CMF results include a trivariate probability 
distribution from the following considerations of three variables, 

(7 = waves, u = upper level, L = lower level), 

1. ® (a, w,) 

2. n(u.w,) 
within each of these figures, six plates (consisting of three pairs) 
are used to show the three possible combinations of the three variables. 
The pairs arise because all results are computed with and without the 
inclusion of random numbers. For example, six plates corresponding to 


the trivariate (1) above (where n(u,w) means JDF of u and w with n as 


CMF) are: 
Withouk Random Numbers With, Randomehunbexs 
1. mn (u,w) 4. n(u,w) 
Zemewin,u) . 5. w(n,u) 
Ss.) u(n,w) 6. u(n,w) 


Phase-amplitude results are presented in two figures, one figure for 
all variables at each level and the wave. Each of these figures include 
the following arrangement of the variables: 


i «WI Top Panel in Figure 18 


Ze UW 

a. WwW 

4. u 

Ss) eal Bottom Panel in Figure 18 
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A. RESULTS FOR 19 AUGUST 

This period was from 1226 to 1244 CST. Measurements in the air were 
made at 1.5 and 4.0 meters above mean water Tevail Higher waves existed 
in this pericd than during the other period considered, 26 September. 
The following descriptions for general conditions were given by Davidson 
(1970): 


"The weather map for 0700 and the hourly wind and wave histories for 
the day appear in Figure 12 (of the present study). A warm front had 
previously passed over southern Lake Michigan earlier in the morning. 
As a result of this frontal passage the winds were southerly and the 
atmospheric stratification was stable. The temperature difference 
between the 4.0 meter level and the surface was 1.2 C to 1.5 C. From 
mid-morning to late afternoon, there were scattered clouds. Wind (at 
the 16 meter level) and wave histories (H1/3 am, hourly averages) 
indicate that the waves began to develop under the influence of the 
steady south wind following frontal passage. Even though the wind 
increased during the day, the wave field appears to have obtained 
equilibrium by mid-morning (0900) and remained so throughout the 
observation pericd. The wind direction and wave propagation direction 
were observed to coincide. 


ae ee we 


Figure 13 shows wind (2 and 8 meters), temperature (degrees Kelvin), 
and wave (Hj, /3) conditions for the observational period. During the 
period there was a slight increase in wind speed but the wave field 
appears to have decreased during the last half of the period. The 
dynamics of the matched layer should have had no influence on the 
results because the critical level was effectively at infinity." 

Variance and covariance spectra for this period, obtained by Davidson 
(1970), appear in Figure 14. Velocity variance and covariance spectra 
for the 1.5 meter and 4.0 meter levels appear on sides A and B respectively. 
The wave spectrum, along with the bandpass filter (dashed line), appears 
on both sides in order to relate the wave spectrum peak to extrema in the 
velocity spectra. 

Primary features in these spectra are extrema in both variance and 


covariance spectra near 0.2 Hz, the wave spectrum peak. These occur for 


both levels. Energy maxima appear near this frequency in the variance 
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FIGURE 12. General Conditions for 19 August 1968 
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FIGURE 13. General Conditions for 1226 to 1244 CST, 19 August 1968 
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FIGURE 14. Velocity Variance and Covariance Results and Wave Spectra 
(With Bandpass Filter Superimposed) for 19 August 1968, 
With 1.5 Meter Level on the Left and 4.0 Mecer Level on 
the Right 


a2 





spectra of both components at both levels near or centered on this 
frequency. Covariance spectrum extrema for 1.5 meters indicate an upward 
transfer of momentum, ¢ | (n)>0, and that for 4.0 meters indicate relatively 
larger negative values for ead (n), resulting in, perhaps, enhanced 
downward transfer of momentum. 

Because these results were obtained at a fixed levei above mean water 
level, energy concentrations at the frequency of the wave spectrum peak 
could reflect simple bending of streamlines, as depicted in Figure 15A, 
in a coordinate system moving at the speed of the surface wave. In this 
realization, observed fluctuations in u and w would occur as shown in 
Figure 15B and correspond to those observed in potential flow. Referring 
to Figure 15B, u and w are in quadrature, 90° out of phase, and would not 
contribute to the shear-induced momentum flux. 

In the foilowing paragraphs, properties of the wave-induced moticn in 
the airflow will be examined in order to determine why contributions to 
the momentum transfer were observed and why opposite contributions were 
observed at two levels along a plumb line. Comparable relations for a 
linear process, shown in Figure 15, suggest that the observed spectral 
results are due to non-linear properties of the wave-induced fluctuation 
within the shear flow. Dynamics of the critical level were not a factor 
in the results for this period because the critical level was well above 
the level of measurements, i.e., the base speed of 0.2Hz wave is equal to 


7 msec”! and the windspeed at 4.0 meters was 5 m sec7!. 
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FIGURE 15. Streamlines for Potential Flow; (A) Simple Bending of 
Streamlines, (B) Relations Between u and w in Potential 
Flow 
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JDF/CMF results for u, w, and the waves (n) appear in Figures 16 and 
17 for the 1.5 and 4.0 meter levels respectively. A general comparison 
between plates A in these two figures reveals a significant difference 
between the u-w JDF contours at the two levels. U-w JDF contours at 1.5 
meters, Figure 16, are close to being concentric circles while those for 
4.0 meters, Figure 17, are asymmetrical with a major axis 135° 
counterclockwise from the positive u-axis. 

Although results appearing in plate A of these figures represent, 
statistically, the probability of joint occurrences of fluctuations of 
various magnitudes of u and w, further discussion of the results will be 
within the context of physical processes. Therefore, these results will 
be interpreted from the point of view that the respective quadrants, 
counterclockwise from the positive u-axis, represent the following 
processes: 

1. Quadrant I - upward transfer of positive momentum; 

2. Quadrant II - upward transfer of deficit momentum; 

3. Quadrant III ~- downward transfer of deficit momentum; and 

4, Quadrant IV - downward transfer of positive momentum. 

Herein, negative fluctuations in u represent deficit momentum and 
positive fluctuations represent positive momentum relative to the mean 
momentum at the observation level. 

General patterns occurring in plate A agree with those physical pro- 
cesses interpreted from the cospectral results. The 1.5 meter cospectrum, 
with the band defined by the bandpass filter, indicated that, in the mean, 
there was a reduction in the net downward transfer of momentum. The reason 
for this (Figure 16, plate A) is probably because there is a decrease in 


downward transfer of positive momentum (Quadrant IV) anc an increase in the 
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FIGURE 16. U-W JDF/CMF Results for 1.5 Meters for 19 August 1968 
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FIGURE 17. U-W JDF/CMF Results for 4.0 Meters for 19 August 1968 





downward transfer of deficit momentum (Quadrant III). The cospectrum 
(Figure 14) for the 4.0 meter level represented enhanced net downward 
transfer of momentum. This is complemented by the JDF results in which 
the probability contours are asymmetrical, elliptical, and represent 
Significantly greater probability for upward transfer of deficit momentum 
(Quadrant II) and downward transfer of positive momentum (Quadrant IV). 

For the 1.5 meter level slight asymmetry of low probability contours 
are indicative of greater than normal occurrences of larger negative 
fluctuations in u at this level. These large negative fluctuations in 
u do not appear to have a preferred association with either positive or 
negative fluctuations in w. The latter interpretation follows from the 
fact that the probability contours have a flat edge appearance at negative 
u values. 

The role of the waves in producing the above general patterns and 
deviations can be examined from CMF contours for the wave heights in 
plates A of Figures 16 and 17. The general CMF patterns agree with those 
expected for measurements at a fixed level within bending streamlines, 
potential flow. Observed fluctuations in u and w with respect to the 
waves for potential flow appear in Figure 15. 

For potential flow, the CMF patterns for the waves would be maximum 
positive wave values associated with negative u and zero w. Maximum 
negative wave values would be associated with positive u and zerow. A 
sufficient deviation from the pattern is the bimodal, or split, maxima 
and, to a lesser extent, minima in the wave CMF distribution. There are 
sufficient points contributing to the split maxima to maintain this 


pattern even with the addition of random numbers (Figure 16B). 
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The statistical significance of the bimodal extrema can also be 
evaluated by examining the JDF's for u and the waves and w and the waves 
which appear in plates C and E of Figure 16. The asynmetry expected in 
potential flow, due to negative correlation between u and n, does not 
sar in plate C of this figure. There exists a very slight asymmetry 
along the 135° axis but it is not as large, for example, as that between 
uw appearing in plate A of Figure 16. In fact, the JDF for u and n are 
nearly as concentric as those for w and n, which in potential flow should 
be concentric because the two variables are 90° out of phase. 

Probable causes for the deviations in both JDF and CMF results, from 
that predicted by potential flow theory, will be examined from phase- 
amplitude results. Phase-amplitude relations for the variables u, w, T, 


“uw, and wI with respect to the wave appear in Figures 18 and 19 for the 


Aan . a ss per ee ; j 
4.0 meter levels respectively. As previously indicated each 


(Qu 


2.5 an 
figure has three panels wnich represent: 

1. small waves (A) - waves less than 1.50; 

2. large waves (B) - waves greater than 1.50; and 

3. all waves (C). 

In Figure 18 the wave appears to be a near perfect sine function but 
the traces for u and w are distorted. Distortion in u and w are most 
evident within the negative fluctuation of each. In general, at both 
levels, extrema in w occur over the node being positive over the rising 
part of the wave and negative over the sinking part. In general, extrema 
in u are 180° out of phase with the wave being negative over the crest 


and positive over the trough. 
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FIGURE 18. Phase-Amplitude Results for 1.5 Meters for 19 August 1968 
(A) Small Waves, (B) Large Waves, (C) All Waves 
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FIGURE 19. Phase-Amplitude Results for 4.0 Meters for 19 August 1968 
(A) Small Waves, (B) Large Waves, (C) All Waves 
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A significant deviation from pure potential flow is the distortion 
in the u phase~-amplitude for the 1.5 meter level. Negative fluctuations 
in u, at 1.5 meters, extend further behind the wave than predicted by 
potential flow theory. In fact, the maximm wees value of u occurs 
behind the crest. 

Distortions which appear insignificant within the phase-amplitude 
depiction for u and w are emphasized in the -uw phase-amplitude results 
which will now be considered. Figure 20 is included in order to compare 
the observed results with those predicted by a linear expansion, i.e., 
Fourier series, of the fluctuations. Phase-amplitudes represented in 
Figure 20 were obtained with sine waves using phase angles, with respect 
to the wave, obtained from the spectral analysis. The resulting symmetry 
along the wave is evident in this figure. Also, there is a mean downward 
momentizn transport (-uwele) when the phasc relation between u and wis 
greater than 90°. 

In comparison to the linear predictions shown in Figure 20, -uw phase- 
amplitudes obtained in this study were asymmetrical. The -uw phase- 
amplitude results, for large waves, panels B in Figures 18 and 19, reveal 
that the possible processes (upward and downward transfer of positive and 
negative momentum) are not as predicted by results from a linear description. 
In fact, downward transfer of positive momentum, denoted as point A, is 
much smaller than any of the other processes. 

In summary, these results, which correspond to a case where the waves 
were in a state of equilibrium or slight decay, indicate that, as a first 
approximation, the airflow adjacent to the waves conforms to the irregular 
mobile wave surface. This was particularly true for the 1.5 meter level 


Where JDF contours were nearly concentric. However, more detailed analyses 
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of JDF/CMF results and also analyses of phase-amplitude results revealed 
Significant distortions, from the potential flow prediction, in the u 
component. This was not the case for the w component which exhibited, 
at both levels, phase relations predicted by potential flow. 

Phase~amplitude results for -uw appeared to show most clearly the 
processes associated with a wave-induced fluctuation in the airflow. At 
both levels upward and downward transfer of deficit momentum appeared 
to be the dominant process. Therefore, assuming that random or shear- 
induced properties are averaged out during phase~-amplitude computations, 
these results indicate that the waves influence did not enhance the 
downward transfer of positive momentum. 

A probable cause for the distortion in the u component, at both levels, 
appears in Figure 21. It is suggested that a high pressure region, 
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centered slightiv behind the crest an 
wave, produced the observed decelerations, and hence distortions, from 

the potential flow predictions in the u component. Such a pressure 
pattern coincides with the general condition of a decaying wave field. 
Because the pressure extrema are moving, it is apparent that this influence 
on the wind field is not the same as that by a stationary object. These 
interpretations suggest that, for this case, the drag coefficient should 
approximate the smooth flow case with some adjustment for the observed 


pressure influence. 


B. RESULTS FOR 26 SEPTEMBER 

This ie was from 1355 to 1413 CST and the measurements were at 
1.5 and 4.0 meters. Although the waves were not as high during this 
observation period as they were on 19 August, they increased slightly 


during the period. The wave spectrum was observed to have changed shave 
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during the day. This change was a shift of the spectral peak to a 
lower frequency. The weather for the day (abridged from Davidson, 1970) 
is shown in Figure 22. The wind during the observaticn period was 
relatively constant and the hydrostatic conditions were unstable with 
the water temperature approximately 5 C higher than the air temperature 
at four meters. Wind and wave histories are shown in Figure 23. 

Variance and covariance spectra for this period obtained by Davidson, 
appear in Figure 24 in which the velocity spectra for the 1.5 meter and 
4.0 meter levels appear on sides A and B respectively. Again the wave 
spectrum peak, with the bandpass filter (dashed line), appears on both 
Sides. | 

Similar to 19 August, the primary feature in these spectra are 
extrema in both variance and covariance spectra, at both levels, near 
the wave spectrum peak approximately 0.26 Hz. Covariance spectra indi- 
cate enhanced downward transfer, by Tt) <9, of momentum for both levels 
contrary to the 19 August spectra which showed an upward transfer of 
momentum at the 1.5 meter level at 0.2 Hz. 

JDF/CMF results for u, w, and n appear in Figures 25 and 26 for the 
1.5 meter and 4.0 meter levels respectively. A general camparison of 
the two levels (plate A) reveals that they are similar. That is, the 
u-w JDF in both have a major axis of asymmetry approximately 135° counter- 
Clockwise from the positive u-axis This general pattern agrees with the 
physical processes interpreted from the cospectral results which suggest 


a net downward transfer of monentum. 
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-W JDF/CMF Results for 4.0 Meters for 26 Sept 


PIGURE 26. U 





It would be difficult to differentiate the JDF results from those 
for a "normal" shear flow. In a shear flow, u-w JDF probability contours 
would also be asymmetrical, elliptical, and would represent significantly 
greater probability for upward transfer of deficit momentum (Quadrant IT) 
and downward transfer of positive momentum (Quadrant IV). An apparent 
deviation is that Quadrant II and IV are not symmetric. Upward transfer 
of deficit momentum (Quadrant II) is the dominant of the two. Wave QF 
features and phase-amplitude relations will be the primary acts in 
identifying the influence of the waves. At the 1.5 meter level, joint 
occurrences between large negative fluctuations in u and positive fluctu- 
ations in w are strongly associated with the wave crest. This vareeion 
is not as evident at the 4.0 meter level. The CMF patterns at both 
levels are, of. course, not the same as those predicted by potential flow 
theory (described in the 19 August discussicn). Similar te 19 August, 
panels C and E do not show sionificant asymmetries between either u or w 
and the waves. However, the lower level, 1.5 meters, does have the higher 
correlations. 

The manner of the deviations from potential flow is quite apparent in 
the phase-amplitude results, Figures 27 and 28. Similar to 19 August, 
the wave appears to be a nearly perfect sine function. At both levels, u 
appears to have the same general relation as in potential flow, negative 
over the crest and positive over the trough. This relationship is almost 
perfect at 1.5 meters. The 4.0 meter level is distorted during the 
positive cycle with the maximum u shifted forward of the trough and has 
a flat appearance. Although, the phase of extrema in u, for 4.0 meters, 
is discernible in these results, the amplitude is very small and, in fact, 


barely above a level which would be viewed as significant. 
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FIGURE 27. Phase-Amplitude Results for 1.5 Meters for 26 September 1968 
(A) Small Waves, (B) Large Waves, (C) All Waves 


(B) 


60 








(A) (B) Ke) 


FIGURE 28. Phase-Amplitude Results for 4.0 Meters for 26 September 1968 
(A) Small Waves, (B) Large Waves, (C) All Waves 
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The primary phase-amplitude result is the shift in the w maximum. 
Whereas in potential flow, and on 19 August, the w maximm is ahead of 
the wave crest (zero over the crest), in this observation period w 
maximum is directly over the wave crest. 

The deviations are again emphasized in the -uw trace and Figure 29 
is included for comparative purposes. Figure 29 depicts the product of 
two sine waves whose phase angle, with respect to the wave, are the 
same aS u and w and which were obtained from spectral analysis. It can 
be seen by comparing Figure 28 and 29 that the momentum transfer processes, 
as depicted by JDF/CMF results, differ significantly from the results, 
depicted by a linear description. The dominate process (large scale 
waves (B) in Figures 27 and 28) appears to be upward transfer of deficit 
momentum and the only other process of any significance appears to be 
downward transfer of positive momentum. Downward transfer of positive 
Momentum is, however, much less than the upward transfer of deficit 
momentum. 

In summary, these results, representing a case in which the dynamics 

£ the critical level and wave generation processes could be a factor, 
reveal significant differences from those observed in the other period 
examined. U-W JDF results for both levels indicated that upward transfer 
of deficit momentum and downward transfer of positive momentum were the 
most significant processes. These same results indicated that the former 
was the most dominant of these two processes. 

Phase-amplitude results indicated that the asynmetrical, shear flow 
type, pattern of the u-w JDF were due to a near 180° phase difference 
between u and w. This phase difference existed because u maintained out- 


of phase relations with the wave extrema but the maxim i in w shifted, 
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FIGURE 29. Phase-Amplitude Results for Sine Waves Utilizing Spectral 
Analysis Phase Angles; (A) 1.5 Meters, (B) 4.0 Meters for 
26 September 1968 








relative to the 19 August description, to a position over the crest. 
The shift in phase of the w component is predicted in the linear theory 
and is due to the dynamics at the critical level. 

The prominence, again, of deficit momentum transfer was evident in 
the -uw phase relation. A physical interpretation for this, as was 
possible for 19 August, has not been achieved. However, it could be 
related to the slight distortions in both the u and w phase-amplitude 
traces. It is known that distortions of this type are associated, in a 
linear expansion, with super-position of different harmonics. It is, 
therefore, suggested that in this case the distortions were due to a 
coupling between the wave-induced fluctuations and a fluctuation of a 
higher frequency. The latter could be due to buoyant processes which 


existed over this relatively warm undulating wave surface. 
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V. CONCLUSIONS 


In this study, turbulence data obtained over waves were examined to 
obtain samé understanding of properties and influence of wave-induced 
fluctuations on the turbulent regime. Statistical procedures were used 
which enabled interpretations to be made on non-linear properties of 
the wave-induced fluctuations and on separate modes of momentum transfer. 
Two data periods, representing decaying and generating wave conditions, 
were examined and the results revealed significant difference in these two 
periods. The usefulness of the selected statistical procedures is 
supported by the following conclusions from the results of each of these 
two periods. 

For a case when the near surface wind speed is much less than the 
speed of the predominant surface wave, the critical level nearly at 
infinity, these results indicate: 

1. The primary feature of the flow is one of simple streamline 
bending in which the airflow conforms to the wave surface and this can 
be observed as high as 4.0 meters above mean water level. 

2. There is a slight deceleration in the airflow due to a pressure 
Maximum located near the wave crest. Because this pressure maximum is 
traveling at the speed of the wave it does not represent an obstacle to 
airflow as a stationary object would. 

3. An observed increase in the negative correlation between hori- 
zontal and vertical components, which is normally considered to represent 
enhanced downward transport of positive momentum, was due to enhanced 


upward transfer of deficit momentum. 
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4. If these results are representative of conditions when the waves 
are in equilibrium with the wind field, an appropriate drag coefficient 
would be one reflecting smooth flow conditions. This is based on the 
fact that the airflow simply follows the wave. 

For a case where the near surface wind speed and the phase speed of 
the predominant surface waves are about the same, the critical level near 
the level of the observations, the results indicate that: 

l. There is an increased negative correlation between u and w that 
is due to a 90° phase shift of the wave-induced w component. This shift 
is predicted by a linear theory which is based on the dynamics at the 
critical level. | 

2. There appears to be a distortion, from potential flow predictions, 
in the u component and this could be due to wave-induced buoyant motion. 

Future research on these data and other data obtained over natural 
water waves should include provisions for the statistical procedure 
employed in this study. Wave following sensors (for exammle Boston, 
et.al., 1969, Davis, 1969) should be used in the measurements. This 
would remove the necessary comparisons with potential flow predictions. 
Relations between temperature fluctuations and the waves should be 


interpreted. 
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APPENDIX A: JDF/CMF COMPUTER PROGRAM 
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PENA Do LOMIPRRALURE JDF/CME ANALYSES 
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FIGURE 31. JDF/CMF Results; 
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FIGURE 32. JDF/CMF Results; Uy 5 and Uy 0 
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FIGURE 34. JDF/CMF Results; Ty 5 and Ty 0 with Waves for 19 August 1968. 
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FIGURE 35. JDF/CMF Results; W 


1.5 and T. 5 with Waves for 26 September 1968. 
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FIGURE 36. JDF/CMF Results; W and T 


4.0 4.0 with Waves fcr 26 September 1965. 
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FIGURE 37. JDF/CMF Results; and U with Waves for 26 September 1°68. 
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FIGURE 38. JDF/CMF Results; Wis and W, ) with Waves for 26 September 1968. 
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